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Abstract. 

Surveys of galaxy clusters provide a promising method of testing 
models of structure formation in the universe. Within the context of 
our standard structure formation scenario, surveys provide measurements 
of the geometry of the universe and the nature of the dark energy and 
dark matter. Cluster catalogues will be constructed using some combina¬ 
tion of X-ray, optical/near-IR, and mm or cm-wave observations. These 
catalogues will be used to study the cluster redshift and mass distribu¬ 
tions along with the correlations of the cluster spatial distribution. These 
measurements probe the volume-redshift relation, the power spectrum of 
density fluctuations and the evolution of galaxy cluster abundance. All 
are sensitive to the amount of dark matter Dm, the amount of dark en¬ 
ergy De, the equation of state of the dark energy w{z) and any other 
parameter, which affects the expansion history of the universe. 


1. Introduction 

Over the last few years, cosmological constraints from Type la SNe (Schmidt 
et al. 1998; Perlmutter et al. 1999), cluster baryon fractions (White et al. 
1993a; David et al. 1995; White & Fabian 1995; Buries &: Tytler 1998; Mohr 
et al. 1999; Arnaud & Evrard 1999), the cosmic microwave background (CMB) 
anisotropy (Hanany et al. 2000; Jaffe et al 2000; Lange et al. 2001) and other 
complementary measures (Bahcall et al. 1999, and references therein) have 
pointed toward a dark energy dominated universe (Da ~ §)) with a significant 
dark matter component (D^ ~ |) and a trace of baryonic matter. The recent 
detections of the 2nd and 3rd acoustic peaks in the CMB anisotropy (Halverson 
et al. 2001; Netterfield et al. 2001; Pryke et al. 2001) lend additional support to 
these conclusions and bring several important questions into sharp focus. At the 
dawn of this new era of precision cosmology, the important questions concern 
the very nature of the dark matter (collisionless or self-interacting) and the 
characteristics of the dark energy (which we can parametrize by the equation of 
state parameter w, where the pressure p = wp). 

Recent theoretical and experimental developments make future cosmolog¬ 
ical studies that utilize galaxy clusters extremely promising. One particularly 
promising approach is the use of galaxy cluster surveys, which enable one to mea¬ 
sure the cluster redshift distribution and the correlations in the cluster spatial 
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distribution. Surveys are now being carried out using cluster X-ray emission, 
the near-IR/optical light from cluster galaxies, the distorted morphologies and 
alignment of background galaxies, and the effect that hot electrons within clus¬ 
ters have on the cosmic microwave background (the so-called Sunyaev-Zel’dovich 
effect or SZE; Sunyaev & Zel’dovich 1972). To use these surveys to full effect in 
cosmology studies, we must first test the standard model of structure formation. 
In addition, we must sharpen our understanding of the nature and evolution of 
galaxy cluster internal structure and the relationships between cluster observ¬ 
ables (i.e. SZE decrement. X-ray emission, galaxy light) and the cluster halo 
mass. 

In these proceedings we describe a fundamental test of the hierarchical 
structure formation model, and then we examine in some detail the cosmological 
dependences of the cluster redshift distribution. We end by highlighting some of 
the challenges that currently exist in using cluster surveys to precisely constrain 
cosmological quantities like the equation of state of the dark energy. 


2. Structure Formation Constraints from High-z Cluster Surveys 

Because of the nature of the power spectrum of density fluctuations, we expect 
that structure formation proceeded hierarchically from small to ever larger scales 
(i.e. Peebles 1993 and references therein). Low mass galaxy clusters ( r\j lO^^M©) 
are expected to first emerge at redshifts of z = 2 to z = 3 within the currently 
favored model. Higher mass clusters (~ lO^^M©) appear later at lower redshifts. 
An appealingly powerful test of structure formation would be to probe the clus¬ 
ter population with sufficient sensitivity to detect the first emerging low mass 
systems. High sensitivity SZE surveys are particularly well suited for studies of 
the high redshift galaxy cluster population, because of the redshift independence 
of the decrement AT: 


AT 

Tcmb 



dlneksTe, 


( 1 ) 


where Tcmb is the cosmic microwave background (CMB) temperature, aT is the 
Thomson cross section, me is the electron rest mass, c is the speed of light, 
/cb is the Boltzmann constant and Ue and Te are the electron number density 
and temperature. In other words, if one has a galaxy cluster described by a 
particular distribution of Ue and Tg, the magnitude of the SZE distortion of the 
CMB along a line of sight passing through the cluster would be independent of 
the cluster redshift. This together with our expectation for how cluster structure 
evolves with redshift, makes SZE instruments capable of detecting clusters of a 
particular mass no matter what that cluster’s redshift (Holder et al. 2000). 
This is a particularly power approach for determining the redshifts when galaxy 
clusters first emerged. 

In contrast to this SZE behavior, cluster X-ray emission (along with any 
emission) suffers from (I -|- z) cosmological dimming. The X-ray surface bright¬ 
ness Ix is 


1 /^e 

dvr (1 -|- z)^ hH 


dlnlA{Te), 


( 2 ) 
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Figure 1. Interferometric SZE and ROSAT X-ray (inset) images of 
three galaxy clusters of comparable mass at redshifts z = 0.2, 0.5 and 
0.8. SZE contours are 75 ^K, and the X-ray color scale is the same 
in each cluster. Note that while the SZE signal remains comparable 
at increasing redshift, the X-ray surface brightness suffers cosmologi¬ 
cal dimming. This characteristic makes SZE observations particularly 
well suited to studies of high redshift clusters, (figure courtesy J.E. 
Carlstrom and J.J. Mohr) 


where Uempge = P, "mp is the proton rest mass, p is the intracluster medium 
mass density, and A is the temperature dependent X-ray emission coefficient 
describing bremsstrahlung and line emission. This strikingly different behavior 
of X-ray emission and the SZE is qualitatively illustrated in Figure 1 by the panel 
of interferometric SZE observations (with X-ray image insets) of three clusters 
of comparable mass at redshifts z = 0.2, 0.5 and 0.8. The SZE contours and 
X-ray color scales are the same for all three clusters. Although the cluster SZE 
signal is similar at all redshifts, the X-ray emission dims rapidly, as expected. 

High sensitivity interferometric SZE surveys carried out with a new gener¬ 
ation of SZE optimized interferometers will soon carry out the fundamental test 
of hierarchical structure formation described above. Three such instruments, 
the SZ-Array, AMiBA and AMI (all described elsewhere in this volume), are all 
funded and currently in various stages of construction. 


3. Cosmological Constraints from the Cluster Redshift Distribution 

Within the context of the standard structure formation scenario, it is possible 
to use cluster surveys to measure cosmological parameters. The abundance of 
galaxy clusters and its redshift evolution have been recognized as sensitive probes 
of the normalization of the power spectrum and the mean matter density in the 
nearby universe (White et al. 1993b; Viana &: Fiddle 1999). The parameter 
degeneracy between ag, the rms amplitude of mass fluctuations in the universe 
hltered on an 8h~^ Mpc scale, and the matter density parameter 11^ can be 
broken by extending cluster surveys to higher redshift (Bahcall et al. 1997). 
This particular probe is highly complementary to the CMB anisotropy, because 
it probes the era of structure formation- when dark energy becomes dominant- 
as opposed to the era of recombination. 
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Figure 2. The cluster redshift 
distribution of an interferometric 
SZE survey similar to the one 
planned with the SZA. All curves 
are normalized to produce the ob¬ 
served local abundance of mas¬ 
sive clusters. The redshift distri¬ 
bution is sensitive to cosmological 
parameters. Note the inset, which 
shows the limiting mass of clus¬ 
ters, which result in 5a detections 
in mock SZE observations. 


The observed cluster redshift distribution in a survey (see Figure 2) is the 
comoving volume per unit redshift and solid angle dVjdz dQ times the comoving 
density of clusters ricom with masses above the survey detection limit Mum- 
written as 


dN 
dz dQ 


dVr., 


dz dQ 


'^com — 


H{z, 


■d\{z){l + z) 


'Mu, 


^ dn 
dM 

^(z) dM 


( 3 ) 


where dn/dM is the cluster mass function, H{z) is the Hubble parameter as a 
function of redshift and d^ is the angular diameter distance. The cosmological 
sensitivity comes from the three basic elements: 

• Volume: the volume per unit solid angle and redshift depends sensitively 
on cosmological parameters (i.e. higher Qa or lower increases the vol¬ 
ume per solid angle). Figure 3 (left) is a plot of the comoving volume 
element {dV/dz/dM) versus redshift for three cosmological models. Note 
the rapid increase in the volume element at modest redshift, which is re¬ 
sponsible for the rapid rise in the cluster redshift distribution in Figure 2. 
At higher redshift the comoving volume element flattens out and eventually 
turns over. 

The cosmological sensitivity of the distance-redshift and volume-redshift re¬ 
lation derives essentially from the expansion history of the universe E{z), 
where H{z) = HqE(z), where EIq is the Hubble parameter and the pa¬ 
rameter E{z) describes its evolution. Within our cosmological framework, 
the expansion history of the universe simply depends on the nature and 
amount of the constituents that make up the universe. That is, E‘^{z) = 
(1 + -|- (1 — ~ ^ e ) (1 T ^) T (1 + z )^^ ^ \ 

• Abundance: the number density of clusters at a given redshift depends 
sensitively on the growth rate of density perturbations. This growth rate is 
highly sensitive to cosmology (i.e. higher speeds the growth of density 
perturbations so that clusters “disappear” more quickly as we probe to 
higher redshift). Figure 3 (right) is a plot of the comoving abundance of 
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clusters above a fixed mass, where the abundance is normalized to reproduce 
the observed local abundance of massive clusters. Note that abundance 
differences increase dramatically with redshift and are responsible for the 
high redshift {z > 1) behavior of the cluster redshift distribution (Figure 2). 

As stated above, the cosmological sensitivity of the abundance evolution 
appears to derive from the growth rate of density perturbations. Within 
the linear regime, the differential equation that describes growth depends, 
again, on the expansion history of the universe E{z). The rapid evolution 
of the abundance is due to an exponential dependence of abundance upon 
the amplitude of density fluctuations on the galaxy cluster scale (Press & 
Schechter 1974; Jenkins et al 2001). 




Figure 3. The comoving volume element (left) and cluster abun¬ 
dance above a fixed mass (right) in three different cosmological mod¬ 
els. The abundances are normalized to produce the observed local 
abundance of massive clusters. Differences in the cluster redshift dis¬ 
tribution are dominated by volume at low redshift and by abundance 
at high redshift. 


• Mass limit: the mass of a cluster, which is just luminous enough to ap¬ 
pear above the detection threshold, typically depends on the luminosity or 
angular diameter distance as well as the evolution of cluster structure- both 
are sensitive to cosmological parameters. The survey yield and redshift dis¬ 
tributions are both sensitive to the limiting mass, as indicated in Figure 4. 
Figure 4 shows the cluster redshift distribution in a fiducial cosmology for a 
limiting mass of M = 2 x 10 ^^Mq, and for limiting masses 10% above and 
below this value. 

The challenging aspects of using galaxy clusters to constrain cosmology, 
aside from building the instruments to carry out the surveys, include an under¬ 
standing of how cluster abundance evolves within a variety of cosmologies and 
how to relate cluster observables like X-ray emission, SZE distortion, galaxy light 
and weak leasing shear to halo mass. These relations are required for all red- 
shifts. Theoretical studies of structure formation suggest that the mass function 
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Figure 4. The cluster redshift 
distribution within a fiducial cos¬ 
mological model for a mass limit 
of M = 2 X IO^^Mq and for mass 
limits 10% above and below this 
value. The mass sensitivity of the 
survey yields and redshift distri¬ 
butions means that accurate cos¬ 
mological constraints require un¬ 
biased estimators of cluster halo 
mass. 


0 12 3 4 

Redshift 

dn/dM{z) is well behaved, and may be described by a “universal” form when 
suitably parametrized (Jenkins et al 2001; White 2001). Further study is clearly 
required. Observational studies of galaxy cluster scaling relations suggest regu¬ 
larity in the cluster population similar to the regularity in the elliptical galaxy 
population (Mohr &: Evrard 1997; Mohr, Mathiesen &: Evrard 1999; Horner, 
Mushotzky &: Scharf 1999). Hydrodynamical simulations of cluster formation 
suggest that scaling relations between cluster observables and halo mass evolve 
in a simple way (Evrard et al. 1996; Bryan &: Norman 1998), even in the pres¬ 
ence of some early preheating (Bialek, Evrard & Mohr 2001). Much more study 
using higher resolution simulations that incorporate additional physics is clearly 
required here to enable more accurate, unbiased estimators of cluster mass. 

A recent study by Diego et al (2001) suggests that a joint analysis of the 
cluster redshift distribution and the observed scaling relations (all available from 
the same survey data) can allow one to solve for the evolving scaling relation 
and cosmological parameters simultaneously (see also Verde, Haiman & Spergel 
2001). More complete studies of the degeneracies between the evolution of clus¬ 
ter scaling relations and cosmological parameters in the analysis of cluster sur¬ 
veys is ongoing. 



4. Precision Cosmology with Galaxy Cluster Surveys? 

Recently, Haiman, Mohr &: Holder (2001) emphasized that large cluster surveys 
extending to high redshift can in principle provide precision measurements of any 
cosmological parameter, which affects the expansion history of the universe (i.e. 
Dm; and w). This is emphasized in Figure 5, which shows the Icr, 2 ct and Su 
joint constraints on and the equation of state parameter w of the dark energy 
for an X-ray cluster survey which yields ~ 10^ clusters with measured emission 
weighted mean temperatures (and therefore virial mass estimates). Only flat 
models (Dm+DA = 1) are considered and the fiducial model Dm = 0.3, Da = 0.7, 
constant w = —1 and h = 0.65 is adopted. Note the Dm — w degeneracy. Also 
shown (dashed line) is the Dm — w degeneracy for CMB anisotropy and SNe 
la distance measurements. The CMB degeneracy assumes that the angular 
scale of the first peak (at fixed h = 0.65) is known to an accuracy of 1%, 
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Figure 5. We plot es¬ 
timated confidence regions 
for joint constraints on 
and the equation of 
state parameter w of the 
dark energy from an X- 
ray survey proposed in a 
NASA Small Explorer com¬ 
petition. All cosmological 
models are flat, and a fidu¬ 
cial model of Qm = 0.3, 
Oa = 0.7 and w = —1 
is assumed. Also shown 
are the confidence regions 
which correspond to a 1% 

-1 - 0.9 - 0.8 - 0.7 - 0.6 ^ r 

^ measurement of the angu¬ 

lar scale of the first acoustic peak in the CMB anisotropy and a 1% 
measurement of the luminosity distance to z = 1. This figure illustrates 
the power of a cluster survey, which yields ~ 10^ clusters with measured 
temperatures. In addition, it shows that the parameter degeneracy in 
the cluster constraints is roughly orthogonal to the parameter degen¬ 
eracy from the CMB or SNe la measurements, (hgure from Haiman, 
Mohr & Holder 2001) 



whereas the SNe la degeneracy assumes that the luminosity distance to z = 1 is 
known to 1%. This figure indicates that an X-ray survey yielding 10^ clusters 
has comparable constraining power to 1% CMB or SNe la measurements. In 
addition, the roughly orthogonal degeneracy between the cluster constraints 
and those from the CMB and SNe emphasizes the complementarity of these 
independent constraints on cosmological parameters. 

The requirements for such precision are (1) a large cluster sample extending 
to intermediate or high redshift and (2) cluster mass estimators that are unbiased 
at the ~5% level. Assuming these two requirements can be satished, cluster 
surveys have as much potential to reveal the nature and amount of the dark 
energy in our universe as either high redshift type la supernovae observations or 
observations of the anisotropy of the cosmic microwave background. 

In fact, with accurate mass estimators, the cluster redshift distribution is far 
more cosmologically informative than simple distance measurements. Figure 6 
contains a plot, which compares the cosmological sensitivity of SNe la distance 
estimates to cluster survey constraints, assuming distances and cluster masses 
are both accurately estimated from the data. The figure plots the ratio of the 
cluster redshift distribution (solid lines) and the luminosity distance (dashed 
lines) as a function of redshift for a few cosmological models. The denominator 
in each ratio is the quantity from the fiducial model, taken to be Qm = 0.3 
and Ha = 0.7 in this example. The degree to which the ratio deviates from 1.0 
provides an indication of the sensitivity to differences in the two models. Except 
for a narrow window around z = 1, cluster redshift distribution contains more 
cosmological information than do luminosity distances. Given the discussion 
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Figure 6. The ratio of the clus¬ 
ter redshift distribution (solid 
line) and the luminosity distance 
(dashed line) as a function of red¬ 
shift for three cosmological mod¬ 
els. The fiducial model is Om = 
0.3 and = 0.7, and the redshift 
distribution and distances from 
this model appears in the denom¬ 
inator of all ratios. Aside from 
a narrow window around z = 1, 
cluster redshift distributions are 
far more sensitive to cosmologi¬ 
cal parameters than are luminos¬ 
ity distances. 


in Section 3 above, this is easy to understand. Surveys probe the volume- 
redshift relation, which scales as the square of the distance. Abundance evolution 
depends exponentially on the growth rate of density perturbations. In addition, 
limiting masses depend on luminosity or angular diameter distances. Naturally, 
these different dependencies can interfere constructively or destructively. 

Of course, cluster surveys will only achieve high precision if cluster masses 
can be accurately estimated, on the average, from observables. That is, precision 
requires that systematic biases in mass estimators be small. This is similar to 
the case with type la supernovae, in that SN la distances are only accurate to 
the extent that the SNe themselves are standard candles. The high potential of 
both approaches has led people to invest significant effort in better understand¬ 
ing possible sources of systematics. Currently, our theoretical understanding of 
the formation and evolution of clusters is less developed than our understanding 
of the dynamics of density perturbations (well within the linear regime) at and 
before the epoch of recombination; however, it seems to me that our under¬ 
standing of structure formation has progressed well beyond our understanding 
of why SNe la’s form a one parameter family of standard candles that have not 
evolved since before the universe was one quarter its present age. This theoreti¬ 
cal heritage in structure formation is an important resource as we move toward 
interpreting ongoing and planned surveys. 

Precision cosmology with clusters requires large cluster ensembles extend¬ 
ing over large (10^-10^ deg^) solid angles. X-ray and SZE surveys of this sort 
will not necessarily have the optical/near-IR data available on every system to 
estimate a sufficiently accurate redshift. Nevertheless, X-ray and SZE surveys 
are attractive, because high signal to noise detections are possible, projection 
effects are minimized, and observations indicate that tight scaling relations in¬ 
volving X-ray observables exist. The price of not having redshifts can be severe, 
as shown in Eigure 7. This figure shows the confidence regions in Qm and w 
corresponding to a 4000 deg^ SZE observations indicate that tight scaling rela¬ 
tions involving X-ray observables exist, survey carried out from the South Pole. 
Confidence regions include marginalization over as, and only flat models are 
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considered. The contours correspond to constraints using only the total number 
of detected clusters, whereas the solid region denotes the constraints in the case 
that cluster redshifts are available. As was made clear in Section 3, the cluster 
redshift distribution is cosmologically rich. 

Possibilities for redshift followup include large solid angle, multiband optical 
or near-IR photometric surveys and direct spectroscopic followup of member 
galaxies. It may even be possible to obtain rough redshift estimates directly from 
the SZE or X-ray data, but further work is required to explore the feasibility 
of this approach (Diego et al, in preparation). Photometric redshift estimates 
of multiple galaxies within each cluster should be sufficiently accurate to allow 
precision cosmology. The required redshift precision is set less by the scale 
of change in the theoretical cluster redshift distributions (which are smoothly 
varying- see Figure 2) than by the fact that accurately inferring cluster masses 
from measured fluxes requires redshifts. 



Figure 7. Confidence regions on 
^1m and mean w from a bolomet- 
ric SZE survey carried out over 
4000 deg^ with a proposed South 
Pole 8 m telescope. This survey 
yields ~20,000 clusters extending 
to redshifts z = 2. This plot un¬ 
derscores the importance of hav¬ 
ing redshifts for detected clusters; 
the solid confidence regions corre¬ 
spond to constraints from the clus¬ 
ter redshift distribution, whereas 
the contours correspond to confi¬ 
dence regions from the total num¬ 
ber of detected clusters. 


The Sloan Digital Sky Survey is one such multiband, photometric dataset, 
and it will be extremely useful in estimating redshifts for clusters at z < 0.6 
(i.e. relatively low redshift samples like the high mass Planck Surveyor cluster 
sample). However, for reasonably deep X-ray surveys and high sensitivity SZE 
surveys, only a small fraction of the sample will lie at these low redshifts, and 
so deeper, multiband followup will be required. The bad news is that the effort 
required to carry out deep, multiband surveys in the optical and near-IR is 
comparable to the effort required to execute the initial SZE or X-ray survey. 
The good news is that there are several projects being designed independently of 
planned and proposed SZE and X-ray surveys that will provide the required data. 
These survey projects include PRIME— a NASA Small Explorer Class proposal 
in Phase A study that will survey one quarter of the sky in the l-3//m range, 
VISTA— a 4 m class telescope with a wide held near-IR and (eventually) optical 
camera that will carry out surveys in the southern hemisphere, the VST— an 
SDSS-like survey telescope with a large optical camera operating in the southern 
hemisphere, and the Large Synoptic Survey Telescope— a 6-8 m class telescope 
with wide held of view to carry out frequent, repeated imaging of large portions 
of the sky. In addition, there are extremely useful large held of view CCD 



10 


Mohr 


cameras available at KPNO/CTIO and on the CFHT. With these projects and 
others together with the exciting science possible with cluster surveys, it is only 
a question of time until very large cluster catalogs can be derived from large 
solid angle, multifrequency surveys. 


5. Discussion 

This contribution contains a description of two ways of using cluster surveys 
to learn about structure formation and cosmology; (1) the context free test 
of hierarchichal structure formation using SZE cluster surveys, which are suffi¬ 
ciently sensitive to detect low mass clusters no matter what their redshift, and 
(2) the use of cluster redshift distributions within the context of our standard 
model for structure formation to determine the quantity and nature of dark 
matter and dark energy in the universe. It’s important to emphasize that there 
is additional information that comes with a cluster survey. This information 
allows one to study the cluster mass function dn/dM as a function of redshift, 
likely improving the constraints derived from integrals over the mass function 
(i.e. equation 3). In addition, surveys (perhaps with some targeted followup) 
enable one to study cluster scaling relations such as the X-ray, optical or SZE 
luminosity-temperature or luminosity-mass relations; a combined study of scal¬ 
ing relations and the redshift distribution may well allow one to solve for the 
scaling relation evolution and cosmological parameters simultaneously (Diego et 
al 2001). 

One can also study the spatial correlations among clusters to infer proper¬ 
ties of the underlying power spectrum of dark matter density fluctuations. With 
good halo mass estimates like those required to use the cluster redshift distri¬ 
bution to full effect, it should be possible to use the cluster power spectrum 
constraints to improve limits on the neutrino mass density. Even in the absence 
of accurate halo mass estimates, it should be possible to use large surveys (in 
volume and number) to measure the scale of the break in the transfer function 
for the evolution of density perturbations. Recently, Cooray et al (2001) have 
emphasized that the physical scale of the break in the transfer function, which 
is the horizon scale at matter-radiation equality, depends on the matter density 
and CMB temperature. The matter density is measured to high precision with 
CMB anisotropy observations such as those with MAP and Planck. Therefore, 
the break in the transfer function is a standard rod, whose scale is independent of 
redshift and is calibrated to high accuracy with CMB data. Thus, measurements 
of the cluster correlation function within redshift shells returns the angular di¬ 
ameter distance as a function of redshift, much like the SNe la but with a strong 
physical basis for the lack of evolution in the standard rod. This approach is 
very complementary to the cluster redshift distribution approach, and it hinges 
less on extracting unbiased estimates of cluster masses from cluster observables 
like the X-ray or SZE luminosity. 
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